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SUMMARY 


In previous work the authors have emphasized, from funnel -cloud-length 
interpretation, that the exceptionally severe tornado is characterized by 
peak swirl speed relative to the axis of rotation on the order of 90 m/s, 
with 100 m/s being close to an upper bound. Further, they have emphasized 
that the thermohydrodynamic achievement of the pressure deficit (from ambient 
pressure) necessary to sustain such swirls requires that a particular struc- 
ture exist. Specifically, a dry, compressional ly heated, relatively slow 
downdraft of air with tropopause-level properties must lie within an annulus 
of rapidly swirling, originally low-level air ascending on a near-moist- 
adiabatic locus of thermodynamic states. This so-called two-cell structure 
is analogous to the eye/ eye wall structure that marks evolution of a tropical 
storm to a hurricane, albeit for the tornado the horizontal spatial scales 
are vastly reduced. The two-cell structure furnishes an observable possibly 
accessible to a passively instrumented, geosynchronous meteorological satellite 
with resolution refined for mesoscale phenomena, for early detection of a 
severe tornado. Accordingly, the low-level turnaround region, in which the 
surface inflow layer separates to become a free ascending layer is 
examined quantitatively; this is the first step in establishing properties 
of the two-cell structure of the severe tornado. From previously established 
properties of the surface inflow layer under the high-speed portion of the 
potential -vortex portion of the tornado (just prior to separation), the 
turnaround flow is satisfactorily modeled as inviscid, such that turbulent 
diffusion is relatively unimportant. Preliminary results indicate that swirl 
"overshoot," i.e., swirl speeds in the turnaround region in excess of the maximum 
achieved in the potential vortex, is modest. 


INTRODUCTION 


In previous work (ref. 1), the authors discussed indicators of the 
evolution of a rotating thunderstorm (also, referred to as a tornado 
cyclone and a mesocyclone) to a storm of even greater severity. The 
relatively long-lived rotating thunderstorm (ref, 2) is characterized 
by a horizontal scale of 10 km and a vertical scale of 10 km, and by 
organized swirl with radial profile of Rankine-vortex type; that is, 
relative to the axis of rotation, the time-averaged azimuthal component 
of velocity increases monotonically with (cylindrical) radial distance 
from the axis to a peak value, and then decreases monotonically with 
further distance from the axis of rotation. The rotating thunderstorm, 
which tends to translate to the right of the mean ambient wind, has peak 
swirl of roughly 20 m/s; it is a vortex of one-cell type : low-level 

inflow, central ascent and upper-level outflow. About one-quarter to one- 
half of these storms evolve to a stage of greater severity in which a 
central, compressional ly heated nonrotating central downdraft is inserted, 
such that the swirling updraft is displaced from the axis of rotation. 

Such tornadogenesis bears some analogy to the insertion of a relatively 
dry eye within an almost moist-adiabatic eye wall, the evolution that marks 
the intensification of a tropical storm to a hurricane (or typhoon) 

(e.g., ref. 3); of course the temporal and horizontal scales are greatly 
reduced for the midlatitude local -severe-storm phenomenon, relative to 
the tropical-cyclone phenomenon. Still, in both cases, the result is that 
dry nonrotating air whose th ermodynamic state is roughly t^^ the. Ambient 
tropopause may descend at modest rate through t he ^epAh^ pf^tJiA JtropAAP,^^^^ » 
and the resulting increased low-altitude pressure deficit of the storm center 
from ambient may sustain an increased rate of swirling of up to about 
110 m/s. In both cases, the evolution to so-called two-cell structure 
(nonrotating central descent within an annulus of rapidly swirling ascent) 
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may proceed rather rapidly relative to the time scale for intensification 
of the one-cell vortex. In both the tropical -storm and local -severe-storm 
cases, about one-third to one-half of the formidable one-cell vortices 
intensify to the potentially devastating two-cell vortices. 

Briefly, the authors (refs. 1,4) proposed that a necessary and 
sufficient condition for the exceedingly dangerous, exceptional tornado 
is the existence of an "eye" within an "eye wall"* (Figs. 1,2). The 
significance of this proposal is that it suggested the first explicit 
dangerous-tornado-indicating observable possibly accessible to a passively 
instrumented, geosynchronous meteorological satellite with resolution 
refined to mesoscal e-phenomena requirements. The instrumentation of 
such a proposed satellite cannot penetrate the cirrus deck of the thunder- 
storm anvil. The "hole" in the cirrus shield generated by the compressional 
heating of descent may furnish a possible unique characteristic of the 
exceptionally dangerous tornado. Of course, the descent in the "eye" may 
not persist through most of the troposphere, so the tornado may not 
achieve extreme intensity; existence of a nascent "eye^’ does not assure 
either development or persistence of an "eye" extending through the depth 
of the troposphere. Also, early detection of an existing tornado, as 
distinct from anticipation of its onset, is the most that can be gleaned 
from the observable. Still, since the exceptional tornado can persist 
for over an hour and can have path length of 200 km, and since one rotating 


■k 

The suggestion that a compressional ly heated central downdraft extending 
the depth of the tropopause could serve as the mechanism by which a pressure 
deficit, consistent with about 100 m/s peak swirl could be achieved 
thermohydrodynamically, was set forth around 1970 by three research 
groups, independently (refs. 5-9). (Cf. refs. 10-15.) 
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thunderstorm can spawn a sequence of several such tornadoes in succession 
(ref. 16), the observable is one that might be exploited with significant 
benefit by a meteorological satellite, as part of an early-warning system. 

The more common tornado, characterized by peak swirl of 50 m/s and 
by duration of ten minutes and by path length of a few kilometers, is of 
one-cell type. Without an ^'eye'', and with only moist-adiabatic ascent 
to produce a lightening of weight (such that a horizontal pressure differ- 
ential is generated at low altitude to sustain swirling), the more common 
tornado is consistent with a pressure deficit on the order of one-hundredth 
of ground-level atmospheric pressure. Such a tornado is not detectable 
by means of the observable just discussed, the '*hole'* in a cirrus shield. 
However, perhaps 98% of the death and destruction wrought by twisters is 
caused by approximately 2% of the total number of tornadoes (ref, 17), 

Thus, it seems appropriate to concentrate on an observable that characterizes 
the exceptional tornado of two-cell type, and which is consistent with 
a pressure deficit on the order of one-tenth of ground-level atmospheric 
pressure. Indeed, it may be speculated that '^eye'^ insertion tends to be 
of periodic nature, such that a so-called supercell (i.e., a multiple- 
tornado-spawning tornado cyclone) often undergoes several transitions from 
one-cell to two-cell structure during its lifespan, which may be on the order 
of six hours (refs. 1,18). 

Observors in aircraft flying over the cirrus shield of some well- 
developed thunderstorms have noted "hole formation" associated with local 
pronounced descending air flow and re-evaporation of condensed water 
substance by compressional heating (refs. 19,20). On the time scale of 
one-quarter to one-half hour, severe tornadogenesis occurred within the 
thunderstorm observed, on several occasions. It is reiterated that a 
nascent "eye" by no means assures that a two-cell structure fills the 
troposphere, so it is to be expected that, on some occasions, one or more 
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moderate- intensity tornadoes follow cloud-top descent. Of course, one-cell 
tornadoes may arise without any such "cloud collapse." 

Recently, returns from dual pulsed Doppler radars have been inter- 
preted to infer that tornadogenesis in a mesocyclone entails evolution from 
one-cell to two-cell structure (ref. 18). In particular, tornado onset is 
suggested to be coincident with the occurrence of a dry weak- reflectivity 
downdraft, within a swirling updraft with high reflectivity. Thus, the 
thesis proposed in ref. 1 is given further support.* (See also refs. 19-22.) 

Establishing the spatial resolution that would be required for satellite 
detection of the downdraft "eye" seems worthwhile, even if such resolution 
currently is orders-of-magnitude finer than the resolution of existing 
instrumentation. The width of the tornadic "eye wall" in the upper tropo- 
sphere, for a severe two-cell tornado, is an output to be ultimately 
extracted from the model in the work reported below. Also, one would like 
to know what properties characterize the "eye" and "eye wall" in the upper 
troposphere, because sufficient thermal anomaly in the compressional ly 
heated core, over ambient temperatures at the same altitude, may be detect- 
able one day. 

Of course, the most direct procedure would be to concentrate entirely 
on the thermohydrodynamics of the upper- tropospheric flow in a two-cell 
model of a tall narrow intense atmospheric vortex. However, the tornado 


★ 

Inferences from radar returns should be drawn with caution. For example, 
radar returns are received first at midtropospheric altitudes from rotating 
thunderstorms, and only later are returns received from low levels (ref. 2). 
The inference is drawn that swirling develops first in midtroposphere and 
only later develops at lower altitudes. But radar depends on reflection off 
condensed water substance for a return, and so absence of a return does not 
mean, necessarily, that no swirling exists at low altitudes when swirling 
is observed in midtroposphere. In fact, while density stratification of 
the atmosphere modifies the applicability of the classical Taylor-Proudman 
theorem (ref. 23), still one suspects swirl is not confined to midtropo- 
spheric levels only. Since the condensation level descends as rotation 
increases (witness the descent of the tornado funnel cloud (ref. 24)), 
because energy is transferred to dynamic motion at the expense of static 
enthalpy, delayed low-level returns are quite plausible. (Cf. ref. 25.) 
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must be examined as a total flow system, in order to obtain the proper fluxes 
of mass, momentum, and heat into and away from the upper-tropospheric flow. 
Thus, one must model the entire tornado; in fact, the logical procedure is 
to "process" the tornado "throughput" analytically in the same sequence as 
the tornado "processes" the "throughput" in the field. Hence, one must 
commence with the low-level swirling inflow, proceed to the near-central 
swirling updraft, and only then undertake the high-level swirling outflow. 

The authors are grateful for helpful technical discussions with Robert 
Costen of NASA Langley Research Center, Hampton, Virginia; his encouragement 
throughout this work is appreciated. The computer programs were written by 
Phillip Feldman; the graphical presentations were prepared by Asenatha 
McCauley; and the editing and typing of the manuscript were carried out by 
Patricia DeCuir. 
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SYMBOLS 


C the peak value achieved by the dimensionless streamf unction ijj 

in the separated surface inflow layer 

2 4 2 

F* the distribution of the quantity (r*v*) as a function of ip*, m /s 
F F*/r*^ 

2 2 

G* the distribution of the total head as a function of m /s 
G G*/V*^ 

2 

g* acceleration of gravity, m/s 

h* height of the surface inflow layer, m; also, thickness of the 
separated inflow layer, m 

h 

I the peak value of 4^ in the surface inflow layer just prior to 

separation 

m (1 +Xq ) , also, the perturbational value of x at small values 

of e 

N number of points used to resolve the separated-inflow-layer structure 

n* distance normal to a streamline, m 

P (p* - pp/p*V*^ 

F dx/dz 

p* pressure, kPa 

Q q*/V* 
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(u*^ + m/s 


q* 

r* cylindrical radial distance from the axis of symmetry, m 
s* distance along a streamline, m 
u* cylindrical radial velocity component, m/s 
u u*/V* 

V* the peak swirl speed in the potential portion of a vortex, m/s 

V v*/V* (also written as v) 

V* azimuthal (or tangential) velocity component, i.e., swirl, m/s 

V v*/V* (also written as V) 

w* axial velocity component, m/s 

X 1 - [l - 

X Trr*/h* 

X r*/rj (also written as X) 

y z*/h* 

z* axial coordinate, in cylindrical polar coordinates, m 
z zVrJ 

a h 

3 h 

r* r*V*, m^/s 
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A 


incremental distance in X 


^ z*/r* 

0 angle giving the orientation of a streamline to a line parallel to 

the z* = 0 plane; also, azimuthal coordinate 

3 

gas density, kg/m 
4) Q/(l - X) 

X r*/r| (also written as x) 

3 

i|;* streamf unction, m /s 

ip r* 

Subscri pts : 

a ambient ground-level value 

1 time node in a finite-difference grid 

j spatial node in a finite-difference grid 

0 starting value; also, pertaining to the outer edge of the vortex 

1 pertaining to the condition at which v* = V* 

Superscripts: 

perturbational quantity, dependent on X only 

dimension "Jess quantity (also denoted by absence of an asterisk) 

* dimensional quantity 

->■ vector quantity 

unit vector 
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ANTICIPATION OF TORNADOGENESIS 


Before details of the "eye"/"eye wall" portion of a mature severe 
tornado are sought by modeling, some further brief* comments on the use 
of geosynchronous mesoscale-resolution meteorological satellites for 
anticipation of tornadogenesis are set forth. 

Typically, a midtropospheric minimum in the total static energy 
(the sum of static, gravitational, and potential -latent-heat-of-phase- 
change contributions), such that large low-level values are recovered in 
the upper troposphere only, is suggestive of local -severe-storm conditions 
in the midlatitudes. In fair weather, the profile of the total static 
energy increases monotonically with altitude in the midlatitudes. In 
moderately cloud-free areas, the mesoscale satellite seems incomparably 
capable of detecting convectively unstable situations by continuously 
measuring through the depth of the troposphere over broad expanses, by 
rapidly processing the sounding data, and by judiciously concentrating 
on suspect areas. The difficulty is that whether, when, and in what form 
the possible severe weather is to be manifested, today largely eludes 
understanding. Some appreciation of the role of vertical wind shear and 
internal gravity waves as severe-weather initiators has been attained, 
but tornado forecasting remains a remote goal. At present, modelers have 
not evolved stability analyses which indicate successfully the conditions 
under v/hich a one-cell structure of a tornado cyclone becomes unstable to 
small disturbances, and might undergo a transition to a two-cell structure. 

Thus, after but brief survey, one tends to doubt the near-term prospects 
of a mesoscale-resolution satellite as a means of anticipating tornado onset, 
and to concentrate on rapid detection of severe tornadoes. 
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TORNADO MODELING: METHODOLOGY, AND A PROPOSED STRUCTURE 


The modeling undertaken here is not an attempt to simulate the entire 
vortex system from the fundamental conservation laws and boundary conditions; 
others (refs. 26-28) are pursuing this course and we think that a contribu- 
tion aimed at the understanding of the phenomena in particular regions within 
the storm would be of greater value than yet another overall study. 

The methodology for modeling adopted here involves a conceptual sub- 
division of the mature severe vortex into logical parts, in each of which 
different physical phenomena dominate. Only the locally significant 
mechanisms are retained in treating each subdivision, and thus approximate 
methods of higher analysis (supplemented by modest computational effort) 
succeed in displaying parametric dependence explicitly and yielding physical 
insight. A composite solution might then be synthesized from the locally 
valid solutions by requiring appropriate continuity of dependent variables 
and fluxes at interfaces between subdomains. Here the four-part subdivision 
first introduced by Dergarabedian and Fendell (ref. 8), to model a mature, 
exceptionally severe tornado in time average as an axisymmetric vertical 
closed two-cell system extending from ground to tropopause, is adopted; 
there is diffusive flux, but no convective flux, of any quantity across 
any boundary. It is important to note that the incorporation of an "eye" 
is the major innovation of ref. 8 on earlier three-part structure proposed 
by Turner (ref. 29) and Barcilon (ref. 30). 

The subdivisions are (Fig. 2): (I) a potential vortex and throughput 

supply, characterized by cyclostrophic and hydrostatic balances, with 
negligible radial influx; (II) a parabolic surface inflow layer, fed by 
very slow sinking from aloft and serving as a source of low-level swirling 
influx because of frictional forces ("teacup effect"); (III) a turnaround, 
"eye wall", and outflow region, in which air rises on a moist adiabat, the 
air ascending quickly enough to transport the condensed water substance 
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upward; and (IV) a nonrotating, slowly recirculating "eye" of dry air. 

The tornado is a once- through system, in which warm moist air in I descends 
slowly into II, the one region of the vortex in which friction dissipates 
angular momentum; the swirling influx in II is turned into a swirling upflow 
in III, the "eye wall" sloping away from the axis with height under conser- 
vation of angular momentum in a variable-density fluid; a larger pressure 
deficit from ambient (achieved by dry-adiabatic recompression in (IV)) is 
available to sustain large swirling speeds in I, owing to the outward slant 
of the "eye wall." The moisture-exhausted air that has been "processed" 
through the updraft annulus is deposited in the upper- tropospheric portion 
of region III, and does not interact with the "unprocessed" air remaining 
in region I. The "processed" air deposited in the upper portion of III is 
too dry to descend, then recirculate upward on a moist adiabat capable of 
sustaining a horizontal pressure deficit; tornado lifespans are too brief, 
anyway: if the air in the upper portions of III descended too quickly, 

compress! onal heating would destroy the radial pressure gradient between 
edge and axis requisite to sustaining high swirl (and the system would 
disintegrate). Suggestions by Deissler (ref. 31) that recirculation may 
occur seem unlikely. Exhaustion of the throughput of I, the lower total- 
static-enthalpy property of originally midtropospheric air of I, topographic 
disturbances, etc., cause eventual decay of the system. However, explana- 
tions for decay based on simple radial diffusion of angular momentum in 
time are incomplete and unsatisfactory. 

It is now widely accepted that turbulent shear flow is not just a 
random "mess"; rather, spatially refined, time-dependent measurements 
reveal that these flows consist of coherent vortical interactive structures, 
such that the mean and instantaneous properties of the flow field may be 
quite divergent (ref. 32). The observation that a tall thin atmospheric 
vortex may contain transient smal ler-scale vortices (so-called satellite or 
suction vortices) is about a century old (ref. 33). That unimodal distri- 
butions may not describe the character of fluctuations from the mean at 
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low altitudes in a tornado may evoke academic curiosity; but emphasis 
(refs. 19,34) on transient departures from axi symmetry appears premature in 
terms of the current status of whirlwind modeling, and seems digressive in 
terms of engineering-oriented concern with tornadoes (such as meteorological 
satellite requirements, nuclear power plant design, etc.). 

THE LOW-LEVEL INFLOW LAYER 

The principal subject of the current study is the turnaround, updraft, 
and outflow subdivision of the tornado, region III. This subdivision has 
received relatively little attention, but examination of it is required to 
characterize the "eye” radius as a function of altitude. 

Before attention can focus on region III, it is necessary to review 
what is known about the near-surface layer containing swirling inflow. The 
flow state at the exit of region II furnishes the flow state at the entrance 
of region III. The published solutions for the nonlinear parabolic boundary- 
value problem describing flow over a stationary surface under an intense 
vortex are far too many to enumerate, but only recently has lucid apprecia- 
tion of the subtleties of the flow been achieved. 

The outer flow in region I is a potential vortex: during spin-up, under 

convectively induced advection, fluid particles move closer to the axis about 
which they possess angular momentum. As soon as the pressure reduction by 
inflow and spin-up balances the pressure reduction achieved by near-moist- 
adiabatic ascent in the core, no further inflow is possible, and a quasi - 
steady state is achieved. Thus, in our model of the mature vortex, the outer 
boundary condition is approximately a potential vortex with negligible radial 
influx (ref. 8). At least as early as 1969 (refs. 35,36), results pertaining 
to a tornado boundary layer were reported in which, under a potential vortex 
without radial flow, the tangential velocity component monotonically decreases 
to zero from the inviscid-flow peak, as the ground is approached at fixed 
radial station. The radial inflow increases monotonically from zero in the 
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potential vortex to a near-ground peak, and then decreases rapidly to zero. 
Perhaps because the calculation is not carried in close enough to the axis 
of symmetry, the peak for the radial inflow at any fixed radial position 
never achieves values above about one-half the local swirl speed in the 
inviscid flow. In any case, these results furnish limited insight into 
boundary-layer structure; that contribution is made largely by analytic 
work now described (Figs. 3-6). 

Carrier (ref. 37) and Burggraf, Stewartson, and Belcher (ref. 38) 
noted that, under an inviscid outer flow with swirl speed Increasing with 
decreasing radial distance from the axis of rotation, the role of friction 
is confined to an ever-thinning sublayer of the near-ground inflow layer. 
This frictionally controlled sublayer is immediately contiguous to the 
ground and, under the inner high-speed portion of the vortex, is exceedingly 
thin; across this thin near-ground sublayer, the swirl is uniformly neglig- 
ible, and the radial influx decreases from its local peak value (achieved 
at the outer edge of the thin sublayer, and virtually equal to the swirl 
speed of the inviscid outer vortex at the fixed radial position of interest) 
to zero (achieved at the ground) under the action of viscous effects. The 
preponderance of the surface inflow layer, which increases in thickness 
very modestly with decreasing distance from the axis of symmetry in the 
high-swirl portion of the vortex, is mostly inviscidly controlled. More 
specifically, at any height above the very thin near-surface sublayer, the 
conservation of radial momentum is a balance of pressure, centrifugal and 
inertial effects such that, at fixed radial position, the sum of the squares 
of the radial velocity component and of the azimuthal velocity component 
(the swirl) is approximately equal to the square of the local swirl of the 
outer potential vortex. The conservation of angular momentum is a balance 
of inertial and transverse-friction effects, but under the high-speed 
portions of the vortex, the role of friction on the continuing evolution 
of the swirl profile is not large. The conservation of axial momentum 
gives axial invariance of the pressure across the entire inflow layer in 
the conventional manner, and continuity reveals that a very small downward 
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drift feeds the radial influx (appreciable because the area over which the 
downward drift occurs is large). Of course, under the outer, slower portions 
of the potential vortex, i.e. , further from the axis of rotation, the inflow 
layer is thinner, and friction is appreciable across the entire layer; 
however, it is the inner, higher-speed portions of the vortex that are 
relevant to the present objective of establishing the nature of the low- 
level inflow entering subdivision III. The point is that angular momentum 
is conserved along all streamlines with finite swirl emerging from region 
II and entering region III, because the flow across the width of the sur- 
face inflow layer (aside from the near-ground thin frictional sublayer), 
to excellent approximation, is inviscid. under the high-speed portion of 
the outer vortex (only). However, because of the role of friction during 
radial influx toward the tornado core in the outer portions of the surface 
layer II, the angular momentum of each emerging streamline differs. The 
angular momentum monotonically decreases from the peak value of the inviscid 
flow of region I for the outermost inflow-layer streamline, to zero for 
the streamline at the edge of the thin viscous near-wall sublayer. The 
extension of these results to a turbulent inflow layer is given by Carrier 
and Fendell (ref. 39), and, recently, a numerical study by Prahlad and 
Head (ref. 40) displays similar results. 

An important implication of the results just discussed is that, for 
radial positions within the distance at which the potential -vortex model 
holds for the inviscid outer flow, one cannot prescribe the radial profile 
of the swirl. That swirl profile, together with the radial and the axial 
flow components, emerge in the solution of the problem. This observation 
seems particularly pertinent to the two-cel 1 -structure severe-mature-tornado 
case, in which the surface boundary layer separates prior to reaching the 
axis of rotation. Then, the inviscid character of the turnaround, in which 
conservation of angular momentum describes spin-up from a different datum 
for each streamline of the turnaround with significant swirl, may produce 
a radial profile of the azimuthal velocity component at fixed height that 
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is not easily anticipated. Thus, work (ref. 41-47) that prescribes some 
sort of rigid-body-1 ike profile for the outer swirl at radial positions 
within the distance from the axis at which the potential -vortex- type flow 
holds, seems inadequate. 

Because the nature of the low-level flow in region II is important 
to the development given below for region III, it is noteworthy that, in 
addition to the laboratory measurements of refs. 36, 46, and 47, there 
are field markings indicating a strong radial influx, without rotation, 
into some tornadoes at ground level. The evidence is a so-called deposition 
line sometimes left behind in mud, sand, wood debris, etc., after tornado 
passage (refs. 48-50); small objects lie inclined toward the line tracing 
the path of the bottom of the axis, on both sides of the path. The 
objects lie pointing toward the trace of the axis, but deflected to some 
extent in the direction of whirlwind translation. Of course, existence of 
markings depends upon the surface traversed, and markings are not always 
easily interpretable, even if they do exist. Further, solid objects falling 
through a tornado may retain motion imparted by higher-level flow, where 
swirl is appreciable. So markings left by fallen objects may not be repre- 
sentative of the very- lowest- level flow field. 

PRELIMINARY REMARKS ON THE TURNAROUND REGION 

The difficulty in treating the "eye wall" annulus, region III, is 
that one has a free-streamsurface problem at both confining streamsurfaces, 
the streamsurface demarcating the "eye"/"eye wall" interface and the 
streamsurface demarcating the "eye walT'/inviscid-potential-flow interface. 
The location of each boundary must be ascertained in the course of solution 
of the boundary-value problem. 

Here, the flow in the turnaround region is formulated first in 
cyl indri cal -polar coordinates, and, in that coordinate system it is an 
elliptical problem. Alternatively, a thin-layer approximation can be 
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adopted, in which case the problem is reduced to an ordinary differential 
equation by averaging some annular properties over the thickness; the 
result is an estimate of radial displacement of the annular centerline as 
a function of "eye wall" properties. Finally, the "eye wall" problem is 
rendered a tractable initial -value-type problem by reformulation in terms 
of streamf unction and orientation coordinates and by numerical integration 
after finite-differencing. 

The most interesting aspect of the solution is the manner in which 
the erupting boundary-layer streamlines run in toward the axis of symmetry, 
and then move away from the axis, with increasing altitude. In fact, the 
solution tends to be periodic with altitude, such that successive inward 
and outward excursions of the "eye wall" annulus provide the formal solution 
of the formulation. However, it is highly probable that this configuration 
becomes unstable and that the description is valid only locally in the 
turnaround region. 
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A BOUNDARY-VALUE PROBLEM FOR THE TURNAROUND 


An inviscid incompressible steady axi symmetric model of the turnaround 
region in a severe tornado is adopted; the cylindrical polar coordinates 
have their origin at the ground on the axis of rotation (and of symmetry) 
(Fig. 7). 

The conservation of angular momentum and of mass may be written in 
the form 



0) 

7T az* 

(2) 

r*w* = r*h* 

^ ^ TT ar* ’ 

(3) 


where v*, u*, vj* are circumferential, radial, and axial components of velocity, 
r* = r^V*, and r^ is the radius at which the pressure in the potential 
vortex is the same as the pressure in the stagnant core. Super asterisk 
denotes a dimensional quantity; no asterisk, a dimensionless quantity. The 
cylindrical radial coordinate is r*; the axial coordinate, z*; the maximum 
height of the surface inflow layer, h*; the streamfunction for the secondary 
flow (involving velocity components u*, w*, only ) , Tp*. The function F(i|;) 
describes the distribution of angular momentum, which, of course, is con- 
served on each streamline. 


Conservation of radial momentum is implied by (p^ is fluid density) 

(§)-■" 


^ ar* ^ 9z* ar* 


_ r* F(i|;) 

*3 
r* 


(4) 


and Bernoulli's equation is 
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(5) 



+ g*,* 

2r*“^ 



The gravitational contribution is negligible within the region of interest 
here and the pressure across the inlet of the turnaround region is constant, 
so (5) becomes 




2r*^ 2r*^ '^d 


( 6 ) 


where is the pressure at r|. Within the potential vortex and in the 
boundary layer beneath it, 

P* - Pf = (py2)v*^ = (py2)(r*^/r^^), (7) 


where pt is the ground-level ambient pressure (given). 

a 

From a crude treatment of the turnaround, to be published elsewhere, 
(h*/r^) = 0.2 for relevant r|,V*, where typical (given) values, for a 
severe tornado, are r^ = 160 m, V* = 100 m/s. In any case, h* and r^ are 
taken as known. Thus, in addition to its use in (2) and (3), h* is employed 
to nondimensional ize the independent variables: 


X = 7rr*/h* , y = TT2*/h*. 

Subtraction of the radial derivative of (5) from (4) gives 



( 8 ) 


(9) 


where subscripts x and y denote partial differentiation. The swirl profile 
F(ij;) is known from conditions holding at x = x-j = Trr|/h* s Tr/h. It is known 
that the flow in region II is such that F is a monotonically increasing 
function of where F(0) = 0 and F(I) = 1, and where the datum i|;(x-| ,0) = 0 
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1s adopted, and from (2) and (8), 


IT 

^ u(x^,C) dc = I, 


( 10 ) 


where u = u*/V*, and v = v*/V*. Clearly, for F linear in ip, (9) becomes 
Helmholtz's equation, and for F quadratic in ip, (9) remains linear. Such 
simple forms for F(i|;) are not detailed replications of (r*v*)/r* in region 
II, but they are potentially useful approximations. 

The pressure along the streamline that separates from y = 0,at x = x-j 
remains at p^, for the model of a nonrotating "eye" isobaric at altitudes 
at which gravity plays no role. Then, from (6), since F(0) = 0, 

= (|r) on t); = 0. (11) 


The pressure along the streamsurface = I, which passes through the 
circle x = Xp y = tt, should be consistent with the pressure in the potential- 
vortex region; hence 


p(x,y) = Pi + 



( 12 ) 


where 


P* 

1 p*y ■ 


(13) 


From (6), (12), and (13), since F(l) = 1, 

= 0 on = I. 


(14) 
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This condition merely states that there is no contribution from the secondary- 
flow velocity components u*, w* to the pressure field, where region III 
interfaces with region I. 

Equations (9), (11), and (14) constitute the boundary- value problem 
of interest, when supplemented by a relation for F(i|i), obtained from pre- 
viously executed analysis of region II. By (1), (2), (6), and the fact 
that w(x-j,y) = 0, the relation for F(i|;) gives implicitly an expression for 
^'(x-jjy), 0 < y < 7T, with which to initiate the analysis (see below). Finally, 
it is anticipated that ijj(x,y), as given by the boundary-value problem is 
periodic in y; this periodicity reflects the fact that the formulation must 
be revised after completion of about one cycle. More explicitly, it is 
anticipated that the locus of any streamline, say ip = 1/2, as y increases, 
decreases to a minimum value of x, then increases in x to recover its 
initial value x = x-j and to achieve a peak value in x, before decreasing 
back to X = x-j to begin another period. It is to be expected that the 
undulating flow so generated will be unstable somewhere along its trajectory, 
but our present job is to find that flow and the stability problem is 
deferred until that is done. Of particular interest are the amount by 
which V exceeds the value associated with x = x-j , = I , and the values of 

X and \p at which the maximum occurs. 

For explicitness, a rough characterization of the function 4^(x-|,y), 

0 < y < TT, is now presented. In the largely inviscid portion of region II, 

(6) may be rewritten as 

u^(x,y) + v^(x,y) = (x^/x)^ . (15) 

From (1), (2), and (15), one obtains 

(x-| >y ) / P \ I/O 

-u(x^,y) = ^ = {l - [v(x^,y)]‘^} = [1 - F(ip)r^\ (16) 
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where 0 £ F(i/;) ±1 for Thus, ip(xj,y) is given by 

^(x. ,y) 

y = /• ^ d. » 

{ {l -[F(u)]}^''^ 

where ij; ^ I as y -> it (cf. (10)), For F(i|>) = 


= y(l- ^) ’ ^ = f v^Cx^.y) = ^ (i - . 

For F(if) = (*/I)^, 

>p = 2^''^sin(y/2), I = 2^^^^ v^(x^,y) = sin^ (y/2) . 


In fact, 

v(x^,y) = 1 - exp (- ^y) , c = 3 , 

might be more realistic, but is far less tractable. 

In summary, if 

Y(x.y) , G'fr) =-^, 

Then the boundary-value problem is, for 0 £ y 1, 


x(^fx) ^ '^yy + 6'(>i/) = 0 ; 




2 

X 



on 'F = 0, 



on 'F = 1; 


( 17 ) 


(18) 


(19) 


(20) 


(21) 


( 22 ) 


(23a) 
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given, with ip periodic in y. (23b) 

One might consider interchanging the roles of the dependent variable 
\p and the independent variable x, and using perhaps some technique within 
the framework of the method of weighted residuals, to extract the desired 
information from (21) - (23). Because the authors consider numerical 
treatment of (2l) - (23) worth pursuing, the formulation has been set forth 
in some detail. However, alternative procedures are used in this report to 
study C21) - (23). 

TREATMENT OF A SEPARATED LAYER WITHOUT STRUCTURE 

A momentum balance in an axi symmetric separated boundary layer without 
structure is examined. The thin sheet (or "eye wall") is the demarcation 
between any "eye," isobaric at pressure p^, and a potential vortex with 
radially dependent pressure given by (12). At r* = r^, p* in the sheet equals 
that in the "eye" (Fig. 8). 

Henceforth in this section the symbol r*(z*) denotes the outside 
surface of the "eye wall" without structure. At the point A, given by 
[r*(z*)j z*] , the principal radii of curvature are |l + [r*' (z*)]^P'^^/r*"(z*) , 
the radius of curvature in a plane containing the axis of symmetry and the stream- 
line on which point A lies, and r*(z*)|l +[r*‘ (z*)]^P^^, the radius of 
curvature in a plane perpendicular to the streamline on which point A lies. 

Super prime denotes ordinary derivative with respect to the argument of the 
function. The velocity component in a plane containing the axis of symmetry 
is denoted q*, while the velocity component (swirl) in a plane perpendicular 
to the axis of symmetry is denoted v*. Thus, if t is a unit vector in the 
plane containing the axis, and 0 is perpendicular to t and refers to the 
azimuthal component in a cyl indrical-polar-coordinate system, then 

V* = q*t + v*e. (24) 
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The force balance perpendicular to the thin sheet at [r*(z*), z*] 
equates the pressure gradient consistent with a potential vortex outside 
the sheet, to the component of acceleration perpendicular to the sheet. 
There are two contributions to the acceleration, each involving the 
square of a velocity component over an appropriate radius of curvature: 


^ /p*\ _ q*^r*" v*^ 

^ W/ [1 + (r*')2]3/2 ■ ,*[i + (,*.) -jl/2 


( 25 ) 


where h* represents a coordinate running across the sheet. 

Variations in velocity occurring across the sheet are not resolved 
(r* z r*V*): 



dh* = r*^ 


/ 


dh* 


[r*(h*)] 


2 

r* V* 



2 


2 2 

A* 





A* ^ f dh* 

y [r*(h*)]^ 



*^dh* 



B* 



dh* . 


( 27 ) 


The potential- vortex form of v* is used in the definition of A*. 
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From these definitions and from (12), to an accuracy that serves 
current purposes. 


V*2 j 



2 2 


2 

g*y*^j<**" 


2 ! 

{ r*2j 

^ *3 

r* 

1 + (r*')^ 

1/2 " 

1 + (r*')^" 

3/2 


( 28 ) 


The following nondimensionalization is introduced: 


X = — 


— T* 


Z = 




_ 2A* o _ 2B* 

- y,* » ^ y»* 

rj r^ 


(29) 


If A* = (h*/2), B* = (h*/2), seemingly reasonable values, then 


a 




Under (29), (28) becomes (a, 3 specified) 


(30) 


1 - ^ 
3^ 


a 


3x" 


x^(l 


x'2) 


1/2 


(1 + x'^) 


3/2 


X 


II 




(31) 


In this translational ly invariant equation, it is taken that x(z = 0) = 1. 
It is shown below that x is periodic in z. Solution is sought for positive 
and negative z, where the boundary conditions are 
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x(0) = 1 

(32) 

x'(0) = >(', given const. > 0. 

(33) 


The boundary conditions preclude odd or even solution for ^ in Sought 
are , the largest value of x, which occurs where x’ = 0, and x” , the 
smallest value of which also occurs where x' ~ 0. In that x” <1, 
there is "overshoot" of the swirl, and, from (12), there is decrease of 
pressure from the value that holds in the "eye", i.e., in 0 < r* < r*(z*); 
it is reiterated that the magnitude of the swirl overshoot and pressure 
decrease are of particular interest. 

While numerical integration of (31) - (33) is required ultimately, 
some preliminary treatment is helpful. Since (31) is translationally 
invariant in Y, phase-plane analysis is introduced (Fig. 9): 

^ = P(x) => — = (34 

dz” dz dz dx dx 


Thus , 


bpp 

■f?2 


(1 + ^) 


3/2 


a 

—3 —2 

x'^ (1 + P ) 


1/2 -2 
X 


(35) 


The slope is infinite at P = 0 and is zero where 


OL 


-3 -2 1/2 -2 

x3(i + p2) X 


-1 = 0 ; 


(36) 


26 



the intersection of the curve of (36) with the P = 0 ray is given by 


a 



+ 


1 



1 = 0 , 


(37) 


where x* = 1 for a = 0, x* > 1 for a > 0. For 1 » a > 0, 
x^= 1+ (a/2) - (3a^/8) + ... 

Although the following development is not pursued to the extent of 
obtaining results, it may be worth noting that if 

H(x) = (1 + , (38) 

then (35) becomes 


If 


-3H' 



T 



» 


then 



aH 



(39) 


(40) 


(41) 


or 




+ E exp 


(s)- 


(42) 
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where E is a const, of integration. From (38), (40), and the definition 
= (dx/dz), one may write formally 


Thus, 


dz 


z = {l + [H(t)]'2| 


-2>l/2 _ d(,-l/2) 


dz 


dz 


|i + 


( 43 ) 


(44) 


where H(t) is given by (42). 

For the special case a = 0, i.e., no swirl, multiplication of (35) 
by x' yields 


B ^ (x - 1)^ , 6 

(1 + x ' 2 ) 1/2 X ’ 

where the constant of integration has been written as (B/m), with 


(45) 


m 2 (1 + (46) 

for consistency. It is recalled that x^( = dx(0)/dJ) is a given positive 
finite constant. Inspection of (45) reveals that ^'is maximum at x = 1, 
so (1 + x^ ) ' is the maximum value of (1 + x‘ ) , whence the symbol m. 

At 7' = 0, 


— + 
X 


1 + 


1(1 _ 1) + [ 1(1 _ i)l 1/2 [2 + 1 ( 
2 m^ ■ [2 ^ m^J [ 2 ^ 


1 - 


1/2 


m' 


(47) 


For so m 1, x merge to unity; this case involves a vertically 

separating surface inflow layer, and hence no overshoot. For x ^ so 
m X remain bounded; 
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this case involves effectively horizontal inflow of the separating surface 
inflow layer, and leads to the minimum value of x~ for fixed 3. For 
8=0.1, 1.37 and x ' = 0.73; for 8 = 0.2, x'^=1.56 and x'=0.64. 

The plausible range of 3 = 0(0.2), from a crude analysis of the turnaround 
to be published separately; it is recalled that 3 is associated with in-plane 
motion. These results suggest what proves to be a general trend: x 

decreases as x^ and 3 increase. 

Finite values for a indicate finite swirl; increasing a yields larger x 
and smaller overshoot. From results of numerical integration, for finite a , 
approach finite values as ^ fixed 3. For the plausible values 

a = 3= 0,2, for x^ = 2, 5, 10, the corresponding values of x“= 0.77, 

0.72, 0,70. Hence, swirl overshoots in the range of about 10% seem plausible, 
but no more; estimates that swirl in the turnaround exceeds the swirl at 
X = 1, 7 = h (in terms of the definitions of (29)) by about 100% (ref. 26) 
are excessive according to this analysis. Further results are given in 
Table 1 and graphical presentation is given in Figures 10 and 11. 

A MARCHING-TYPE TREATMENT OF THE TURNAROUND WITH STRUCTURE 

A model problem is adopted to study in more detail the local configuration 
within which transition occurs, from radial swirling boundary-layer inflow to 
vertical swirling "eye-wall" upflux, for a steady, largely inviscid vortex. 
Slightly idealized distributions for the radial and circumferential velocity 
components at a particular radius are adopted as input. The model problem 
requires that the streamline which is closest to the axis of symmetry, and 
which constitutes the "eye"/"eye wall" interface, is adjacent to a stagnation 
region of constant pressure, denoted pp and determined in a more complete model 
by hydrostatic considerations in the vertical; here, p* is taken as given. 

It is recalled that at radius r^, at which the surface-boundary-layer 
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"inlet flow" to the turnaround is prescribed, the pressure associated 
with the outer vortical flow of region I is also P^- ' Again, of particular 
interest is the largest amount by which the speed in the transition region 
exceeds that at r* = r|. It is reiterated that since the flow emerging from 
the boundary layer is largely inviscid, the analysis for the turnaround is also 
to be inviscid. 

Conservation of mass and momentum are written in unconventional 
form because it is convenient to adopt as independent variables (coordinates) 
the streamf unction for the velocity in the (r*,z*) plane, and the angle 9 
giving the orientation in that plane of the streamline relative to a line 
in the plane parallel to z* = 0 (Fig. 12). 

The streamfunction \p* is defined by 


r*q* 


dip* 

9n* ’ 


(49) 


where r* again denotes the conventional cylindrical radial coordinate at 
which the velocity is to be calculated, q* is the velocity in the (r*,z*) 
plane, and n* is the distance normal to the streamline. The streamfunction 
of (49) is identical, of course, with the streamfunction for which 
(see (2) - (3)) 


r*u* 


diP'* 


(50) 


r*w* 


dip* 

dr* ’ 


(51) 


where u* and w* are the radial and axial velocity components, respectively, 
and z* is the axial coordinate. 

The circumferential velocity v* obeys the angular-momentum-conservation 
requirements (see (1)) 

r*^v*2 = F*(ip*) , (52) 
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with determined when the velocity distribution of. v* at r* = rj 

is prescribed. The principal curvatures of a surface = const, are 
(80/8s*) (in the (r*,z*) plane, where is distance along the streamline) 
and sin 6/r*. 


The conservation of momentum along a streamline is equivalent, in the 
inviscid theory, to the statement that the total head is a function of i(j* 
only, i.e. , 


2 2 p* - 2 • 


(53) 


From comparison of (53) and (5), q*^ = u*^ + w*^; like F*(iJj*), G*(i()*) is 
implied in the specification of q* and v* at r* = r*, and = const, 
is a case of particular interest. 

The conservation of momentum in the direction normal to a surface 
(p* = const, is implied by (cf. (25)). 


3 /p* \ _ v*^ sin 0 +2 90 . (54) 

9n* \p 3 / ^ 9s* 

This equation must be rewritten for the coordinate system mentioned earlier. 
It is noted that 


9n* 9n* dip* ^ dip* ^ 


and also that 


90 1 _ cos 0 

9s* " (9s*/ 90) ” (9r*/90) * 


Thus, 


9 / ^\ ^ 1 / v*^ sin 0 , g*^ cos 

3/\Pd/ *'*'>* V ''* ' (Sr*/3d) J ■ 


(55) 


(56) 


(57) 
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Int57), and subsequent equations, it is implied that 3( is taken 

with 0 held fixed, and vice versa. 

Finally, note is taken of the geometric identity 


so that 



6* fixed 


sin 0, 


3r* _ sin 0 
9i|j* r*q* 


(58) 


(59) 


Equations (52), (53), (57), and (59) constitute a set of four equations 
from which the four unknowns v*, q*, p*, and r* can be calculated, where 
each unknown is a function of ip* and 0. 

The appropriate initial conditions are 


p* = p| , r* = r| , F*(i|;*) given, G*(i|;*) given on 0 = 0^ = 0 , 


0 < ij;* < 

~ ^ ~ ^max 


(60) 


A particularly interesting case is F* = (const.) ip*, G* 

with ip*^^, calculable, 
max 

The boundary conditions are G*(i|;* ) = V*^/2 : 

I- 1 1 let A. J 


(const. ) 

^ ^ ^max 


= 0,6 ) = 


Pf , 


(61) 


- Pi = 


Pd 

2 


* I- 


v*^(ip* , 0 ) - V*^ 

^^max’ ^ 


pjr*" r 


(62) 


'*"^(lp* ,0) 

'^max’ ' 


where r* e r|V*, with rJ,V* (and hence p^) given 
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The following nondimensional ization is adopted: 


V = v*/V*. Q = q*/v*. X = r*/r* , C = z*/r* . P = (p* - p*)/pgv*^ , (63) 
F = G = G*/V*^. (64) 


Then, (57), (59), (52), and (53) become, respectively, 

2 

p = V sin 6 Q cos 9 
K - XXg 


Y - sin 0 

ip XQ ’ 

w2 - EM 

V 2 

£ + vl+ p - GM 

2 2 2~ 


(65) 

( 66 ) 

(67) 

( 68 ) 


At r* = r|, the thickness of the layer is denoted h*. Under the nondimen- 
sional ization, (50) becomes, at x = 1, 


II = -u = [g( 4-) - F(i)j)]^^^ , 


(69) 


from (53), with = u*^, p* = p*. Hence, at x = 1, if h = h*/r* , 


, = 

0 |^G(i/^^) - F(i|^^)j 


1/2 






■ / 


max 




0 “ F(^l)J 


1/2 


(70) 


33 



The case of interest, from results for the surface inflow layer, is G(ijj) = 1. 
Here the special case F(i1j) = considered: 


^max di|j^ 





hence. 


^max = 


h 

2 * 


(72) 


In general, for choice of G{i|;), F(i^), (70) gives as a function of h. 

Henceforth, ip is taken as known and is denoted by C. 

^max 

The boundary conditions are, from (61), (62), and (63), 


- 1 

P = 0oni[i = 0; P = 5 — on = C 

2X^ 


(73) 


The initial conditions at 0 = 0 are, for 0 < ip < C, 


P = 0, X = 1. V = Vj<p), Q = [g(i(<) - 


(74) 


For 0 < 0 « 1, solution is sought in the form 


QCii-.e) = Q(,{tp) 

+ 9^/^ q(i|)) + ..., 

(75) 

P(i|j,e) = 

+ 6^''^ p(ip) + .... 

(76) 

v(^,e) = V^M 

j^l - d^''^ r(tp) - ... j , 

(77) 

x('P’S) = 1 

+ 0^*^^ r(4i) + .... 

(78) 
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Substitution of (75) - (78) in (65) - (68) gives, since (67) 
becomes an identity. 


P 


I 



( 79 ) 


r 


I 


'\j, . _ 

0 =?> r(i|j) = const = m ; 


(80) 


Thus, 


q - 


2 ^ 
r + 

0 


P 


0 . 






% 2 
q + m. 


(81) 


(82) 


From the first of (73), p(ijj = 0) = 0; from the second of (82), q(0) = 0 
since V^(i|^ = 0) = F(4( = 0) = 0. Differentiation of the second of (82), 
and substitution from the first of (82), give 



- - Q- + m(V^) ‘ 
m ^0 ^0 


a- 

q 



Qo 


+ m V 


2 

0 


(83) 


(84) 


But at Ip = C, Q = 0, so ~ inspection of (85), 

C 

V? = 2 J Qg(<|)j) , (85) 

0 

since V (C) = 1. For G(i|;), F(i|j) of interest, (84) gives q 0 as 

Oj ^ Oj 

ip C. The quantities r, q, p have now been determined, so (75) - (78) 
may be used to initiate the solution. 


35 



It is helpful to introduce 


Q(4> 


.0) = (l - fy'^%(1',0). 


( 86 ) 


where 4>(C,0) is neither singular nor zero. It is also helpful to 
introduce the following change of one of the independent variables: 


X = 1 - [1 - ■I' = c[l - (1 - x)^]. 


(87) 


Thus, 


^ 2C[1 - (i/C)] 1/2 X) 3X > 


( 88 ) 


Q(x,0) = (1 - X) <f. (x,e). 


(89) 


Equations (65) - (68) are transformed as follows: 


1 „ _ y2 sin 0 ^ (1 - X)2 (t, cos 0 


^x 


<(>X 


XXg 


(90) 


2C ^X 


sin 6 

X<P 


(91) 


F[l|:(X)] 


(92) 


(1 - X)2 *2 w2 , „ . S[>I'(X)] 

2 2 2 


(93) 
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The initial conditions at 


0=0 are, for 0 < X < 1, from (74), 


P = 0, X=l, v= F[. 1 -(X)] = Vg[.^(X)], 

{g [.;>(x)] - F[.MX)]}^/^ 

* = • 

The boundary condition on X = 0 is, from (73), 

P(O,0) = 0. 

If G[ij;(X)] = 1 for = X = 0, then since F[i|j(X)] = 0 
from (93), 

c()(O,0) ^ 1. 


Hence, 


X(O,0) X^(O,0) = 2C sin 0 [x(0,e)] 


X 


The boundary condition on i|; = C is now developed 
of (68) with respect to ip gives, with (67), 


00 + F* (ilj) _ + p _ G' (ip) 

.2 ^3 > 2 * 


^ 2X^ 


X^ 


From (85), 


QQ 




= - (p [(f) + 2(1 - ] 


C' 


(94) 

(95a) 

for = X = 0, 

(95b) 

4C sin 0. (96) 

Differentiation 

(97) 

(98) 
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At ^ C, 


J_ *2 + F'(C) . p _ G'(C) 

2C* 2x' ■ x' 


(99) 


But as ip Ci 


= F(C) sin 9 _ p 

x' ■ x^Q ■ ^ ’ 


( 100 ) 


SO 


-^/(c,e) 




2X 


(101) 


„ cf'[^(x = D] 

(t> (X = 1,0) = — , CG' [il)(X = l)] 

^ / V 


If G' [i|^(X = l)] = 0, 


X (X = 1,0) 


P CF'[i|<(X = l)] 

/(X = l,6) = 


X^(X=1,0) 


(102) 


(103) 


Xx(x = i,e) 


2C^^^ sin 0 
F' [>|)(X = 1)] 


(104) 


The boundary-value problem is now of initial-value, or marching type, 
a significant simplification over (9), (11), and (14). If the solution 
for P(X,0) V(X,0), (j)(X,0), and x(X,0) 'is known at 0 = 0 ^., then upon 
assignment of F[i|i(X)] and G [i|^(X)] , the solution may be found at 
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One might proceed via adoption of an explicit finite-difference treat- 
ment of (90), by means of a four-point grid element with central spatial 
differencing and with special provision for enforcement of the boundary 
conditions (96) and (104) at X = 0,1. One may then obtain cj)(X,6^.^^) from 
(91) by simple differentiation V(X,0^^l) from (92) by simple substitution, 
and P(X,0^.^j) from (93) by simple substitution. The procedure may then be 
repeated to obtain values for the dependent variables as a function of X at 
the next increment in the time-like variable 0.j+2* Unfortunately the 
accuracy of this fully f inite-differenced procedure proves to be insufficient. 

Instead, the method of lines (refs. 51-53) is used, such that only 
the derivatives with respect to x av'e finite-differenced; as a result, coupled 
nonlinear ordinary differential equations, each of the first order, must 
be integrated forward in 0 from specified initial values. Specifically, 

(97) is written, by use of (91) and (92), as 


dXj(9) 

~dT^ 


2C^sin^0(l - X^)^ 

(Xj)^(Xx)j 



(x,)^ 


j = 1.2 


..,N 


(105) 


where Xj(e) = x(Xj,e). with x^^ = 0, 1, and Xj^^ > X j . The X^ are 

chosen, along with N, the number of points used to resolve the structure in 
X for each 0. For the special case G[ip(X)] = 1, from (91)-(93), 




(Xj)^ 


(X)()j ^ 4C^sin^0(l - Xj) 


(Xj) 


(Xj)' 




, 'tc^sin^ed -X.)^ (l 


(106) 


(Xj) (Xx)j 


_(Xx)j] 


Three-part difference formulae linking )^.^j(e), Xj(o)» Xj_i(o) give 

(^ ) . and (xyy^i» thereby couple the differential equations for 
X J XX J 

[dx-(9)/U0] , j = 1, 2, ..., N. From knowledge of x^*(9)j one may extract 

J J 

4>.(0), V.(0), and P.(0) by the procedure alluded to in the previous paragraph. 

J J J 
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The integration in 6 is initiated at 0^ where 0 < 9^ « tt/2; the 
series expansion gives Xj(0-]^). The boundary conditions at and Xj^ are 
enforced in the expressions adopted for ^^XxH’ ^^X^N’ ^^XX^N 

(ref. 53). 

Here the simplest conceivable procedure is adopted as a first trial. 

The parameter N = 3; in particular, 

pl/2 

X^ = 0, X 2 = 1 - X 3 = 1 (107) 


for equal mass flux about X^ for the special case F[i|;(X)]= ijj/C => X(2 - X). Thus, 

( 108 ) 


(Xx>2 


(X3 - X2)(X2 -X^)^ + (X2 - Xi)(X3-X2)^ 


'"XX>2 


(X 2 - Xj)(X3 - X2)(X3 - Xj) 

2 ( x ^ " X 2)(^2 ~ ^^2 ” 


(Xp - X^)(X3 - X2)(X3 - X^) 


(109) 


Also, simple two-point differencing is used to express the boundary conditions 
in terms of X 2 ^ 0 ) • 


[x^(o, 0 )]j( = 4C sin e ^ Xi(s) = l[x 2 (e)]^ - 4 C(X 2 - Xj) sin ( 110 ) 


[x(l.e)] = 2C sin 6 => X 3 ( 0 ) = X 2 ( 6 ) * ^^(X^-X^) sin 6 . 


(Ill) 


The entire problem then degenerates to numerical integration of one nonlinear 
ordinary differential equation for X 2 ^o) 9 = 0^ to 0 = (tt/2). 
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A particular case of interest, as previously noted below (60), is 


G[)j^{X)] = 1 ; 

F[«KX)] = I 2 X(2-X)=^ F'[i|)(X)] = C'^ . 

It follows that, from (74), 

0.(1) ■ (‘-tr ■ 

from (85), 

m = ± C 

From (75), (84), (86), (112), (113), and (114), for small 6, 
Q[(4'/C),e] = [1 - <)>[(K-/C),e] ^ [i - (>|//c)]^/2 


(112) 

(113) 


(114) 


(115) 


+ 0 


1/2 


1 Wc) 

2 c r 

m ^ - 1 

(<l'/C)] 

([1 - (<P/c)]^^^ 

II _ 1 

m L 

^ ■ 2 J 


Jr 


(116) 


where 

(j.[(i()/c) ->• i,e] + 1 + e^'^^(-m) . 

From (76), (82), (112), (113), and (114), for small 6, 

(il^/C) 


(117) 


P[(VC),e] ^ 


\C 


1 - 


(118) 
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The results of this paragraph may be expressed in terms of X 
to furnish initial conditions at small but finite 6 , with which to com- 
mence the numerical integration. 

Integration proceeds from 0 = 0 to 0 = {tt/ 2), where x(O, 0 ) 
achieves its minimum value; the solution may be continued by reflection 
from 0 = (tt/2) to 0 = tt: e.g.. 


X[X, (tt/2) - a] = X[X, (tt/2) + a] (1|9) 

where 0 5 a < (tt/2). The solution for 0 = 0 to 0 = (tt/2) is initiated 
by adopting the negative root in (115), such that X decreases with 
increasing 0 , for fixed value of X. Solution for this half of the total 
"cycle" pertains to the inward excursion of the separated surface inflow 
layer. 

The solution is completed by letting 0 = Tr-a, 0£O£ (^ 1 / 2 ), 
in the boundary-value problem posed by (90) - (96), (112) - (113), The 
dependent variable X achieves its maximum value at X = 1, a = (tt/2); 
the solution may be reflected about a = (tt/2), to complete the solution 
for the "turnaround cycle." Formally, if one identifies a with 0, the 
solution for the outward excursion of the turnaround (discussed in this 
paragraph) may be obtained from the identical formulation that yields 
the solution for the inward excursion of the turnaround (discussed in the 
previous paragraph); the only difference is that one adopts the positive 
root in (115) in the initializing solution when one seeks results for the 
outward excursion. 

Results obtained by the method of lines to the marching- type boundary- 
value problem for the turnaround region with structure are given, for swirl 
profile linear in the normalized streamfunction, and for total head 
invariant with the normalized streamfunction, in Figs. 12-15. A key result 
is the very modest overshoot (^ 10 %) of the swirl during inward excursion 
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of the separated inflow layer. Although the solution warrants considerable 
refinement, the locus of the low-level "eye wall" is now delineated. For- 
mulation to extend the solution for the "eye wall" properties, in a tractable 
manner, to the midtroposphere and upper troposphere is now the next logical 
challenge before modelers of severe tornado structure. 


CONCLUDING REMARKS 

Over about the past decade, the authors have used relatively simple 
analytic methods to study quantitatively a quasisteady, axisymmetric 
model of a mature severe tornado, taken as a vertical closed system. The 
procedure has been to adopt locally appropriate simplification to obtain 
results for subdomains of the vortex; when the results have been obtained, 
a composite solution should be synthesized from the results for the separate 
subdomains by requiring suitable continuity of dependent variable and of 
fluxes at interfaces between subdomains. The methodology was the first, or 
among the very first, to yield the following results on tornadoes, some of 
which are now widely accepted. 

1. The peak swirl speed in a tornado, relative to the axis of 
rotation, never much exceeds 110 m/s; the length of the funnel cloud 
may be used to infer, both readily and quickly, the peak swirl, probably 
as accurately as more tedious techniques, 

2. The pressure reduction from ambient necessary to sustain such 
a large peak swirl speed may be achieved only if the vortex has two-cell 
structure, i.e., if the vortex has a central dry nonrotating "eye" 

of compressionally heated, descending air, situated within an annulus 
("eye wall") of swirling air rising on a near-moist-adiabatic locus of 
thermodynamic states. 
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3. The mature severe tornado consists, then, of a four-part 

structure: a potential vortex with the "throughput" supply of warm 

moist air, in a cyclostrophic balance; a near-ground inflow layer, in 
which friction dissipates angular momentum, and into which fluid sinks at 

a slow rate from above; a turnaround, "eyewall," and outflow layer, in which 
the swirling low-level influx separates from the ground, rises through the 
vertical extent of the troposphere, and moves radially away from the axis 
of rotation in the upper troposphere; and the warm "eye," with modestly 
recirculating flow. 

4. The surface inflow layer under the high-speed portion of the 

vortex consists of a two-part structure: a very thin sublayer, immediately 

contiguous to the ground, in which frictional effects are significant, 

and a much thicker region in which an inviscid description suffices. Under 
a favorable pressure gradient, the thin viscous sublayer becomes thinner 
as the layer moves inward toward separation into an updraft annulus, while 
the inviscid portion grows in thickness but modestly with decreasing 
radial position. In the inviscid portion, the squares of radial 
and azimuthal velocity components sum to the square of the local potential - 
vortex speed, at fixed radial position; hence, the large swirl component 
becomes a large influx component as the ground is approached at fixed 
radial position, until friction enters very close to the ground to enforce 
the no-slip boundary condition and to bring the velocity to zero at the 
ground. 

5. The turnaround region, in which the separating surface inflow 
layer becomes a swirling updraft, is described adequately as an inviscid 
layer, such that turbulent diffusion is not a central contributor to the 
dynamic balance {except in a very thin sublayer immediately contiguous to 
the "eye"). The amount of "overshoot," i.e,, the amount by which the 
swirl in the "corner" flow exceeds the peak value achieved in the potential 
vortex, is modest. 
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The next challenge facing such a modeling approach for application to 
the problem of identifying the occurrence of a mature severe tornado, is 
restoration of hydrostatic effects, such that the turnaround solution initiated 
here is continued into the midtroposphere and upper troposphere. Only 
by so extending the present line of research can properties of the "eye"/ 

"eye wall" interface be delineated. This interface, the unique charac- 
teristic of the severe tornado, furnishes the one observable potentially 
accessible to a passively instrumented, geosynchronous meteorological 
satellite with mesoscale resolution, for the purpose of detection of the 
exceptionally dangerous whirlwind. Detection of two-cell structure in a 
vortex system is far more significant than detection of a Rankine-vortex- 
type, one-cell structure of a mesocyclone; existence of a mesocyclone by 
no means necessarily implies the onset of a severe tornado event, but 
evolution to two-cell structure within a mesocyclone is very highly sug- 
gestive of a severe tornado event. 

These statements perhaps have implications for the most profitable 
direction for tornado research aimed at anticipation of the onset of a 
severe tornado, as opposed to early detection of an already developed, intense 
whirlwind. Specifically, rather than apply linear stability theory to a 
two-cell vortex model in order to establish propensity to multiple- 
vortex (i.e., satellite-vortex) structure, one might better examine the 
conditions under which transition of a one-cell structure to a two-cell 
structure is preferred in the swirling system. Indeed, this particular 
stability problem seems to have been neglected, relative to either the 
mul tipi e- vortex problem just mentioned, or the oft-examined question of 
initial organization of a weak depression that ultimately evolves to 
a mesocyclone. 


In view of (1) the insight and convenience derived from explicit and 
largely closed-form results, (2) the absence of detailed field measurements, 
and (3) the uncertainties of modeling turbulent diffusion in rapidly 
swirling atmospheric vortex, the authors believe that such simplified 
subdivisional modeling of the tornado continues to be the most productive 
path to delineating whirlwind properties and structure. 

TRW Systems and Energy 

Redondo Beach, California 90278 
31 October 1978 
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Table 1. Numerical Results for Extrema of the Turnaround, from Integration 
of the Initial-Value Problem (31)-(33). 


a 

6 

a/B 

x' (0) 

^calc 


^calc 

(")a=0 

0.1 

0.1 

0.0 

2.0 

0.791 

0.789 

1.264 

1.266 

0.0 

0.1 

0.0 

5.0 

0.764 

0.751 

1.326 

1.329 

0.0 

0.1 

0.0 

10.0 

0.742 

0.738 

1.348 

1.352 

0.0 

0.1 

0.1 

100.0 

0.732 

0.727 

1.366 

1.372 

0.0 

0.2 

0.0 

2.0 

0.718 

0.714 

1.392 

1.397 

0.0 

0.2 

0.0 

5.0 

0.671 

0.664 

1.489 

1.497 

0.0 

0.2 

0.0 

10.0 

0.656 

0.647 

1.524 

1.533 

0.0 

0.2 

0.0 

100.0 

0.644 

0.632 

1.553 

1.565 

0.1 

0.1 

1.0 

2.0 

0.820 

0.789 

1.295 

1.266 

0.1 

0.1 

1.0 

5.0 

0.777 

0.751 

1.348 

1.329 

0.1 

0.1 

1.0 

10.0 

0.762 

0.738 

1.366 

1.352 

0.1 

0.1 

1.0 

100.0 

0.750 

0.727 

1.381 

1.372 

0.1 

0.2 

0.5 

2,0 

0.748 

0.714 

1.422 

1.397 

0.1 

0.2 

0.5 

5.0 

0.695 

0.664 

1.510 

1.497 

0.1 

0.2 

0.5 

10.0 

0.677 

0.647 

1.541 

1.533 

0.1 

0.2 

0.5 

100.0 

0.663 

0.632 

1.567 

1.566 

0.1 

0.5 

2.0 

2.0 

0.876 

0.846 

1.212 

1.181 

0.1 

0.5 

2.0 

5.0 

0.841 

0.818 

1.244 

1.223 

0.1 

0.5 

2.0 

10.0 

0.828 

0.808 

1.254 

1.237 

0.1 

0.5 

2.0 

100.0 

0.818 

0.800 

1.264 

1.250 

0.2 

0.1 

2.0 

2.0 

0.845 

0.789 

1.325 

1.266 

0.2 

0.1 

2.0 

5.0 

0.796 

0.751 

1.369 

1.329 

0.2 

0.1 

2.0 

10.0 

0.779 

0.738 

1.384 

1.352 

0.2 

0.1 

2.0 

100.0 

0.765 

0.727 

1.396 

1.372 

0.2 

0.2 

1.0 

2.0 

0.774 

0.714 

1.452 

1.397 

0.2 

0.2 

1.0 

5.0 

0.716 

0.664 

1.530 

1.497 

0.2 

0.2 

1.0 

10.0 

0.695 

0.647 

1.558 

1.533 

0.2 

0.2 

1.0 

100.0 

0.678 

0.632 

1.581 

1.566 

0.2 

0.4 

0.5 

2.0 

0.685 

0.615 

1.652 

1.606 

0.2 

0.4 

0.5 

5.0 

0.617 

0.550 

1.789 

1.772 

0.2 

0.4 

0.5 

10.0 

0.594 

0.528 

1.838 

1.832 

0.2 

0.4 

0.5 

100.0 

0.578 

0.509 

1.877 

1.887 


Notes: Results for (x“) ^ from (40). Also, x for x' = 100 are within 2% 

a=0 0 

of values for x^ for a, 3 studied. 
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THE THREE-PART STRUCTURE OF A MATURE QUASISTEADY 
MODERATE AXISYMMETRIC VORTEX 



igure 1. A schematic, not to scale, of a tornado, with 40-m/s peak 
swirl relative to the axis of symmetry. The radial profile 
of the swirl, at midtropospheric altitude, is like that of a 
Rankine vortex. The peak pressure deficit from ambient, at 
ground, is about 3 kPa. Asterisk superscripts to denote 
dimensional quantities have been omitted. 
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THE FOUR-PART STRUCTURE OF A MATURE QUASISTEADY 
AXISYMMETRIC VORTEX (CLOSED SYSTEM) 



Figure 2. A schematic, not to scale, of the structure of a tornado, with 
100-m/s peak swirl relative to the axis of symmetry. The radial 
profile of the swirl, at midtropospheric altitude, reveals a 
nonrotating central "eye" joined to a nearly potential vortex. 
The pressure deficit from ambient, at ground level, may reach 
about 10 kPa; the outwardly sloped "eye"/"eye wall" interface 
permits fluid in the potential -vortex portion of the tornado 
to achieve swirl speeds consistent with this pressure deficit. 
Asterisk superscripts to denote dimensional quantities have 
been omitted. 
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THE BOUNDARY - VALUE PROBLEM 


CONTINUITY {ru)^ + (rw)^ = 0 

(rV)^ - (rv)^ - (ru)^ r / v i , \ 

RADIAL MOMENTUM = [j^(ru) 2 ] ^ - u(ru)^ 

ANGULAR MOMENTUM ^(rv)^] ^ - u{rv)^ -wirv)^ = 0 


w(ru)2 


BOUNDARY CONDITIONS: z ^ v ^ T/r, u -»■ 0 

z = 0: u = v = w = 0 
r = r^: u = 0, v = r/r^ 


V = + V, 


mol ‘'eddy' ‘'eddy 


0 < z < z^(r) 


5 *{r) 


-r 


ko(r/r)5*{r), 


(l-^)dz; l<2(r/r)6*(r) = K^[z^(r)] 


zi(r) < z < oo 


U2^[r, z^(r)] + y^[r, z-,(r|] 


1/2 


SOLUTION: ru = -rj3(r)g'(7j). rv = rf(r 7 ), w = W(? 7 } 

T7 = z/s(r), s = - r)/(rQ + r)] = ^ - (r/to)^ 

g '2 = 1 _ ^2 ]' + gf' = 0 

f = g = 0 atT 7 = 0; f-^1 as q ^ 


} 


MUST PICK f'(0),77i; W(t?) FOLLOWS 


Figure 3. The third-order parabolic boundary-value problem for time- 

averaged description of the turbulent surface inflow layer under 
a potential vortex over a flat smooth nonrotating boundary. 

An approximate similarity solution holding over the bulk of the 
thickness of the layer, under the high-speed portion of the vortex, 
is developed formally. This solution must be supplemented in a 
thin sublayer, contiguous to the wall, in which the no-slip 
condition on the radial inflow is enforced. Asterisk super- 
scripts to denote dimensional quantities have been omitted to 
simplify notation. 
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THE PHYSICS OF THE INFLOW LAYER 


v{r, z) 


u(r, z)^ + y(r, z)^ + w(r, z) z 


r = r r + zz 
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Figure 4. A schematic of the profiles of the velocity components over the 
entire width of the turbulent boundary layer under the intensely 
swirling portion of a potential vortex. The radial inflow 
achieved near the ground is almost as large as the swirl above the 
inflow layer, at the same radial position. Asterisk superscripts 
to denote dimensional quantities have been omitted. 
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VERBAL INTERPRETATION OF INFLOW LAYER 


FAR FROM THE AXIS OF ROTATION (r^ r < tq) 

RADIAL MOMENTUM : 

FRICTION, INERTIA, PRESSURE GRADIENT BALANCE 
ANGULAR MOMENTUM : 

FRICTION, INERTIA BALANCE 
CONTINUITY : 

SMALL DOWNDRAFT 

CLOSER TO THE AXIS OF ROTATION: (r^ < r ^ Tq) 

RADIAL MOMENTUM : 

INERTIA, PRESSURE GRADIENT BALANCE: 

(VERY CLOSE TO WALL, FRICTION ENFORCES NO-SLIP 
CONDITION) 

ANGULAR MOMENTUM : 

FRICTION, INERTIA BALANCE (FRICTION SMALL) 
CONTINUITY : 

SMALL DOWNDRAFT 


Figure 5. Further summary of the results of the surface inflow layer analysis. 
Asterisk superscripts to denote dimensional quantities have been 
omi tted . 
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Figure 6. Another schematic diagram, not to scale, of the postulated four- 
part structure of an idealized mature severe tornado, again taken 
conveniently as a vertical closed axi symmetric two-cell system 
extending from ground to tropopause. The arrows indicate the 
secondary flow (radial and axial velocity compoennt streamlines). 
The radius of maximum wind (swirl) in I is r^; the outer radius of 
the closed vertical vortex system, r*; the maximum depth of the 
inflow layer, h*. The shaded area suggests where diffusional 
transfer of momentum is significant. 
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Figure 7. Schernatic, not to scale, of an inviscid model of the turnaround 
region of a severe vertical quasi steady axi symmetric tornado. 

The "eye", region IV, is isobaric at pressure p|, over the vertical 
extent of interest here, where p* is also the pressure at r* = rj, 
z* = h* (because the pressure is approximately invariant across 
the surface inflow layer II). The "eye wall", region III, is demar 
cated by two streamsurfaces , i]j*(r*,z*) = const, (the position of 
each to be determined), encompassing the mass efflux from region II 
It is recalled that a potential vortex holds in region I. 
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Figure 8. Schematic of the location of an inviscid "eye wall without 

structure" demarcating the interface between an isobaric non- 
swirling "eye" at pressure pj and a potential vortex. The sheet 
representing the "eye wall" nas thickness h* 0; in the turn- 
around region, its displacement from the axis, as a function of 
height above the ground plane, is denoted r*(z*). 
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Figure 9. ^hase-plane properties of equation (3^), where P = (dx/dz), with 
X the dependent variable (x > 0) and z the independen^t variable. 
Isoclines of zero and infinite slope are noted. At x = 1, the 
distance from the axis of symmetry at which the surface inflow layer 
separates, ^ finite positive slope is adopted. For a > 0, x* > 1; 
for a = 0, X* = 1. The sketched trajectory (solution curve) is 
a limit cycle (closed curve indicative of periodic behavior). The 
periodicity is not of physical interest. 
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Figure 11. The solution to the boundary-value problem posed by equations (32), (33), and (35), 
for a = 0.2, 3 = 0.2, and = 2, given in Figure 10, is further discussed. The 

"pressure" denotes the value of (1 - x"^); the "swirl" denotes the value of 
— 3 —'2 1/2 

ax [1+x ] ' ; the "radial/axial" denotes the value of "swirl" minus "pressure", 
i.e., the left-hand side of (35). 



Figure 12. This schematic of the inward, then outward, excursion of the 

separated inflow layer with structure, indicates the bounding stream 
surfaces = 0 and ip* = C. The in-plane radial-axial flow 
speed is denoted q*; the swirl about the vertical axis, v*. The 
inward excursion need be solved for streamfuncti on- i ncl i nation angle 
0 = 0 to 0 = {tt/2), since the solution for 0 = (tt/2) to 0 = tt 
may then be obtained by reflection. An analogous statement holds 
for the outward excursion. The "eye" is isobar ic at pressure p?, 
but the potential -vortex pressure varies with radial position 
r* as noted. 
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Figure 13. Solution to the boundary-value problem posed by equations (78), 

(80), (105), (107), (112), and (115), for C = 0.1, as obtained 
by the method of lines. The notation is that, for example, 

Xn'(6) = x(X4»0)j where X* is defined in (87), (107). Here the 

U J J 

dimensionless radial position x» defined in (63), is given for the 
inward excursion (m < 0) for 0 = 0 to 0 = (tt/2), and for the outward 
excursion (m > 0) for (it - 6), for 6 = 0 to 0 = ('n'/2). 










p2» m < 0 
Q2 » ffi ^ 0 





Figure 16. The solution for the dimensionless secondary-flow (in-plane, 

radial -axial ) speed Q, defined in (63), for the problem discussed 
in Figure 13. 
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